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ABSTRACT 

Wc develop an effective model for a neutron star with a magnetosphere. 
It takes into account the electromagnetic torques acting on the magnetic 
dipole, the friction forces between the crust and the core, and the gravitational 
corrections. Anomalous electromagnetic torques, usually neglected in a rigid 
star model, play here a crucial role for the alignement of the magnetic dipole. 
The crust-core coupling time implied by the model is consistent with the 
observational data and other theoretical estimations. This model describes the 
main features of the behavior of the magnetic dipole during the life of the star, 
and in particular gives a natural explanation for the n < 3 value of the breaking 
index in a young neutron star. 

Subject headings: stars:neutron - stars:magnetic - stars:rotation - pulsars:general 
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1. Introduction 

Pulsars are identified with rapidly rotating, highly magnetized neutron stars. The 
basic information we receive from this kind of stars is a sequence of electromagnetic pulses 
with a very stable frequency, which is interpreted as directly related to the star rotation. 
The angular velocity Q generally decreases gradually as a result of the torque exerted on 
the star by the radiation reaction. In the vacuum dipole model the star magnetic field is 
assumed to be a magnetic dipole M that forms an angle 6 with the rotation axis, so that 
the star loses energy by electromagnetic radiation due to its rotation ( [Pacini 1967, 196E ; 



Gunn fc Ostriker 1969| ). This turns out to be the main source for the star energy loss. In 



this simple model the evolution of the angular velocity for a neutron star with a momentum 
of inertia / is 

^ = -3-^^!^^' (1) 

or, more generally, Cl = —^Q^. Several characteristics of the dynamics of the /i parameter 
have been observed. After glitches a sudden increment of fi has been noticed, which does 
not completely relax back ( [Link, Epstein fc Baym 1992| ; [Link fc Epstein 1997[ ). There is 



also evidence that for old pulsars, with an age of ^ ~ lO^yr, /i is smaller than for younger 
pulsars ([Ruderman, Zhu fc Chen 1997[) . The first behavior can be interpreted as an increase 
of the external torque after de glitches, whereas the second one suggests a slow decrease of 
the torque with the age of the star. 

Besides this, the breaking index n = ^ has been measured in four young pulsars, 
and it takes values between 1.4 and 2.8. If yU were constant the breaking index would be 
exactly equal to 3, but the measured values are smaller than this canonical value implying 
an increasing /i. A number of factors might affect the breaking index. One of them is the 
presence of mechanisms of loss of energy different from the dipolar electromagnetic radiation 
that could change the f2 exponent in eq. (0). For example, pure multipolar electromagnetic 
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radiation gives n > 5 , gravitational quadrupole radiation n = 5 ( [Manciiester fc laylor 



|1977|) , and early neutrino emission n < (|Alpar fc Ogelman 1990|) . These effects are 



expected to be relatively weak for the pulsars with a measured breaking index, and in 
particular the first two would increase n and thus are unable to explain the observed values 
smaller than 3. Another possible factor that could be considered is the effect of the pulsar 
magnetosphere ( |Goldreich fc Julian 1961] ), but although it would change appreciably the 



angular dependence and the numerical coefficients in eq. (|I]), it leaves the Vt exponent equal 
to 3 ( pood fc Ng 1985| ), and thus does not affect the value of the breaking index provided 
that /i remains constant (see references ( [Michel 1983, 1985|) for an alternative). According 
to this, the most natural explanation seems to require variations of the parameters /, M and 
9, which change the breaking index in the following way for a star without magnetosphere: 

VLii n ( e M i\ 

n = 3 + ^- = 3 + ^ 2 + 2— - - . 2 

n\ tanO M I J ^ ^ 

Given that {n — 3) is of the order of the unity, the parameters responsible for the variation 
should have characteristic evolution times of the order of the pulsar age 'Tq = which is 
of around lOOOyr in the case of the stars of known breaking index. 

With respect to the momentum of inertia, sensible variations during such times due 
to changes in the structure, sphericity, or in the coupling of the internal superfluids to the 
crust, seem to be unable to explain the observed deviations ( [Link fc Epstein 1997|) . 

Significant changes in the magnetic moment also seem to be inhibited by the high crust 
conductivity, which can keep the magnetic fields unchanged and thus fixed to the crust 
along times of order lO^yr or more (jLamb 1991[ ; [Chanmugan 199^ ; Phynney fc Kulkarni 



1994[ ; [Goldreich fc Reisenegger 1992[) and hence the breaking index n is unaffected in 



young pulsars. Furthermore, the decay of the magnetic field reduces the radiation rate and 
hence would lead to n > 3. However, there is no general consensus about the behavior of 
the intensity of the effective magnetic moment and other hypotheses have been proposed. 
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Muslinov and Page ( |1995|) consider magnetic fields trapped under the surface during the 
birth of the star, which could rapidly be freed and would thus increase the total magnetic 
moment (see also (piandford, Applegate fc Hernquist 1983| ; [Chanmugam fc Sang 198^ )). 



Ruderman (|1991a, 1991b, 1991c| ) considers the possibility of crust cracking, reordering due 
to vortex pinning and stresses. If fiux tubes in the core proton superconductor are pinned 
to neutron superfiuid vortices, the latter could drag the magnetic structure to the equator 
as the star spins down ( [Ruderman, Zhu fc (Jhen 1997| ; [Srinivasan et al. 19901 ). In this work 



we adopt the first hypothesis and simply assume a stable external magnetic field rigidly 
fixed to the crust, except perhaps in very short periods at the glitches. 

Under the above assumptions, the only responsible for the anomalous breaking index 
could be the variation of the angle between the magnetic dipole moment M and the rotation 
axis O. In this case eq. (|]) becomes: 

n-3 = 2^-^. (3) 

The angular momentum of the star decreases during its life by radiation, and thus Q is 
negative. Therefore a value n < 3 requires an increase in time of 6, e.g. 6 > 0. In brief, 
for a rigid star in vacuum the observed dynamics is related to the behavior of the 6 angle. 
The observational data suggest that its evolution is controlled by three characteristic times. 
There is a fast dynamics, associated to the reorientation of the magnetic dipole during the 
glitches, with a characteristic time less than or of the order of 50 days, an intermediate 
dynamics that dominates during the early stages of the stars, where the magnetic moment 
slides towards the equator and gives the breaking index values smaller than 3, and finally 
a slow dynamics, where the magnetic moment slides towards the rotation axis of the star, 
which involves characteristic times of the order of tq ~ lO^yr ( Lyne fc Manchester 1988 ). 



The intermediate dynamics competes with the slow one, and the observational evidence 
shows that for young pulsars it dominates, leading to the anomalous values of the breaking 
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index, whereas for old pulsars {tq ~ 10^ years) it is the slow one which dominates, giving 
place to the alignment of the magnetic dipole with the rotation axis. However, if we only 
took into account the radiation reaction of the magnetic dipole on a rigid star, the evolution 



of 9 would perform a slow alignement during the life of the star (|Pandey fc Prasad 1996|) 



with a breaking index n > 3, in contradiction with the features mentioned before. 

In this work we consider a more realistic description given by a simple effective model, 
where the star is considered as constituted by two rigid interacting components in presence 
of a magnetosphere. One of them is the core which contains the bulk of the mass, with a 
momentum of inertia Iq and an angular velocity flo- Its angular momentum is mainly given 
by the vortices of superfiuid neutrons. The other component is the crust, with a momentum 
of inertia Ic <^ lo and an angular velocity fie? and with a dipole magnetic moment M. The 
evolution of the system is governed by three kinds of torques: 

a. - torques acting on the crust due to electromagnetic interactions, 

b. - friction torques between the core and the crust due to the interaction of the neutron 
vortices, the proton flux tubes and the electrons in the star, 

c. -torques due to gravitational effects. 

Although this model could seem a rather simplified approach, it is sufficient for 
obtaining a good qualitative comprehension and can be considered as a first step to a more 
sophisticated description. 

In the following section we present a review of the different effective torques we take 
into account in our model. In Section III we write the system of equations that describes the 
dynamics of the model and we analyze the resulting behavior. Finally, the last two sections 
are devoted to the discussion of the results and their phenomenological implications. 
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2. The effective torques 

In this section we review the torques that act on both components of the star. In 
principle the coordinates of our model are represented by three vectors, M, fig and flc, 
which are the magnetic moment, the core angular velocity and the crust angular velocity 
respectively. It should be remarked that the two components considered in the present 
work are identified with the crust and core, and they are not the ones used to describe the 
long relaxation time after glitches, which is currently identified as an effect of the crustal 
superfluid ( [Anderson fc Itoh 1975| ; [Alpar, Cheng fc Pines 1985| ; [Link fc Epstein 199(j| ). Once 



we have chosen a direction of reference there remains five degrees of freedom, because we 
have assumed the magnetic dipole moment constant in modulus and fixed to the crust. So 
the equation of motion for the magnetic field is simply 

M = ricXM. (4) 

In the following subsections we describe the torques that arise from the different interactions. 



2.1. Electromagnetic torques on the magnetic dipole 

We have a magnetic dipole M which forms an angle 6 with the angular velocity flc oi 
the crust. There are several dipole and quadrupole torque terms that act on the magnetic 
dipole, and as we have assumed that it is bound to the crust, these torques are directly 
applied to the latter. They have been carefully discussed by Good and Ng (1985). The 
effect of the average torques that involve quadrupole terms are strongly suppressed, and 
hence it is enough to consider only the dipolar components. There are two types of dipolar 
torques that act on the crust, the non-anomalous torques of order , which are of the 



same order as the classical spin down torque ( pavis fc Goldstein 1970[ ; [Michel fc Goldwire 



1970[ ), and the anomalous torques that are times larger than the former ones, where R 
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is the radius of the star. In the following we will use an orthogonal system of reference with 
z in the Clc direction, y in the plane (fic, M), and x orthogonal to such a plane, according 
to Fig. 1. The y, and z torques, which cause the loss of energy and the alignment of the 
magnetic dipole in a rigid star, are non anomalous. In this plane the torque is 

Ty^^ IcUy.Mx {M xnc)-IcCOyz^c , (5) 

where M = M~"^M, cUyz = | ^yz and ujyz = | ^ i^yz ■ The Uyz and Uy^ coefficients 
contain information on the magnetosphere contribution to the torque. The first term 
changes the magnitude and direction of fic, whereas the last one alters only its magnitude. 
For a neutron star without magnetosphere we have Uyz = 1 and Uyz = 0, but the presence 
of the magnetosphere may greatly change the value of these coefficients. Its effects are in 
general of comparable size to the ones in the vacuum case, and could depend on the angle 
between M and ftc- Their values are related not only to the currents flowing in the near 
magnetosphere, but also to the ones in distant regions, and so cannot be calculated at 
present. For this reason we will maintain these coefficients of the order of the unity as 
unknown adimensional functions of the angle between the dipole magnetic moment and the 
angular velocity, to be phenomenologically estimated. 

The X component has anomalous and normal terms, but the first ones dominate 
because the normal terms are much smaller than the anomalous ones. The anomalous 
contribution is given by 

= Ic cos 6' QcXM , (6) 

where lUx — |^7~p^ (^TTr) Here the Ux coefficient contains information on the 
magnetosphere, and reduces to Ux — 1 when it is absent. In the case of the anomalous 
torque the magnetospheric effects depend mainly on the near region and can be evaluated 
under reasonable assumptions. For example, if we assume that there is a dominant 
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contribution from the Goldreich- Julian charge density on the closed lines of forces, it leads 
to = -1/4. 

As long as ~ 10^ — lO'^ for most pulsars, lj^ is much greater than Uy^ and Uyz. 
Despite this, in the rigid dipole model the anomalous term does not contribute to the 
alignment of the magnetic dipole with respect to fic or to the spin down of the star. 
However, in the two-component model considered here the situation changes and the 
anomalous term acquires a significant role. 



2.2. Dissipative effects in the superfiuid hydrodynamics 

In the outer-core region of a neutron star there is a mixture of superfiuid neutrons, 
superconducting protons and normal electrons. The effects of viscosity and friction 
due to the scattering of the electrons by the neutron and proton vortices give place to 
effective torques between the core and the crust. The resulting crust-core friction can be 
characterized by a time parameter r/. There are several phenomeno logical estimations on 
the basis of the glitch dynamics. For example, upper bounds are given for the Vela pulsar 
of Tf < 10s for a crust initiated glitch and of rj < 440s for a core initiated glitch ( [Abney 



Epstein fc Olinto 1996| ). These values are deduced from the December 24, 1980 Christmas 



glitch data ( |lVlc (Jullock et al. 1990|) . 



This point has also awakened a great deal of theoretical interest. The interior plasma 
couples the neutron superfiuid due to mixing superfiuid effects ( |Alpar, Langer fc Sauls| 
|1984| ; [Alpar fc Sauls 1988|) , and it is locked to the crust due to Alfven waves or cyclotron 
vortex waves. An extensive analysis of these phenomena has been performed by Mendell 
( |1991a, 1991b| , 1997) on the basis of the Newtonian superfiuid hydrodynamics, generalized 



to include dissipation. The characteristic times related to all these couplings are of the 
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order of the second. The mutual friction torque simplifies to 

Tf = fin,-no), (7) 

when the angular velocities differences are small. When the system is taken out of 
equilibrium such as in a glitch, the time response frequency is given hy ^ = cuf = /j^, 
where Ic and Iq are the moments of inertia of the components, and It = Ic + In is the 
total star moment of inertia. Therefore this torque is relevant for the rapid nucleus spin up 
during glitches. The effective friction could depend on the angular velocity. Non dissipative 
effects give torques of the same kind as the gravitational dragging, which is to be discussed 
in the next subsection. If they are also characterized by a time of the order of the second 
they should be smaller than the gravitational torque for a rapidly rotating star. 



2.3. Gravitational effects 

The neutron star is essentially constituted by superfiuid matter, and hence the observed 
rotation can be achieved only by the presence of vortices. The gravitational fields are strong 
enough to be relevant. They induce a change in the shape of the vortex lines, and also 
affect the density of vortices ( [Casini fc Montemayor 1997| ). The main contribution of these 
effects is a correction on the coefficients of the dissipative torques of the order of 15% with 
respect to the flat space-time values, but they do not introduce new terms. 

An additional term arises from the gravitational dragging, which gives place to a new 
torque if the components of the star have different angular velocities. If we neglect the 
gravitational radiation this interaction is conservative and the corresponding torque is 

Tg = zn,xn,, (8) 

where z ~ IJo-^^ — (■^) ^ 0.1 Jc for small velocities, VtR <ti c, where R is the 
star radius and Rg and Rg are the gyration and Schwarzschild radius of the star. This 
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torque has a significant modulus, and at a first glance it could have sensible effects on the 
dynamics of the magnetic dipole. 

The torques Tf and could be seen as the first terms of the mutual torque in 
an expansion around the point of equal angular velocities, and as such the form of the 
interaction is largely independent of the model. We are only considering terms up to the 
first order in the angle between Oc and Oq. As we will discuss later this is enough for our 
analysis. In the following section we will use all these torques to construct a set of equations 
of motion that defines the dynamics of the star components. 



3. The equations of motion 

The equations of motion for the system described in the preceding section are: 

M = fic X M , (9) 

f z 

Oc = UxCos9Qc>^M+ujyzM X {M x nc)-<^yz^c- — i^c-^o) f2c X fio , (10) 

O; = y- (^^c-^^o) +f a X . (11) 

If / ^ the crust decouples from the core and if / — > oo the two components act as a rigid 
body. In both cases we recover the dynamics of the model of a rigid star with a magnetic 
dipole, with momentum of inertia Ic and It = Ic + lo respectively, which has been already 
discussed in the Introduction. 

The eq. (|^) implies that the magnitude of M is constant. Therefore we have only five 
variables in the system, three angles, a, /3,and 6, defined in Fig. 1, and the two moduli of 
the angular velocities, Vt^ and VLq . For these variables the equations of motion are: 

VLc = — — (f2c — ilo cos a) — iOy-Slc sin^ Q—ujySlc , (12) 

/ 

Qo = — (f^c COS a — fio) , (13) 

J- n 
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^i^^^ = ^^(cos/3-cosacos^)+Wj,, sin^^cos^ + ^fio (f2o. (fieXM)) ,(14) 
(i(cos/3) Z^^(^og^_^og^^osa) + fi,fl-f) (Clo.(Cl,xM)) , (15) 



(it /n fi„ \ /, 



cose / / ^ XX . ( f Qo . f 

(cosp — cosacosfc'j . (lb) 

sin a 

In principle this system of equations would be very difficult to solve, but 
in fact it contains several dynamics with very different time scales, given by 
ujyz ~ ojyz ^ ujx <C min(ci;/ ~ j-,fic), which greatly simplify its treatment as we 
will see now. In the ffist place, from eqs. (0) and ([13|) , if Qc and Qo are very different at a 
given instant, they will attain equilibrium in a relatively short time of the order of tj. The 
last equation tells us that the transient of a is also characterized by Tf, and thus, after a 
time of this order, this variable will acquire a value of the order ^ ^ 1- Hence, after this 
transient we will have fic ~ , « ^ 1 and therefore 6 ^ j3, all of them satisfying a slow 
dynamics with characteristic frequencies of the order of uj^ or ujyz. 

As was already commented, the ffist two equations imply that the moduli of the 
angular velocities of the crust and the core will rapidly reach an equilibrium regime where 
(ic ~ do and ~ — ~ (pyz + ^yz siii^ 0^ ^ , i.e., the difference between the angular 
velocities is quickly suppressed. The common decrease of both components is given by 

ilo^dc^ - J- ioJyz + ^yz SVO? O) Vtc , (17) 

which coincides with the one obtained in the dipole model for a rigid star when there is a 
magnetosphere. This means that the torque component along x will not affect the decrease 
of the angular velocity. This is analogous to the case of a rigid star, where this torque does 
not cause any spin-down. It means that the dissipative effects of the core-crust interactions 
do not change the star energy significantly, and therefore the main mechanism of energy loss 
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is still due to the electromagnetic radiation. In fact, the ratio between the power lost by 
friction and radiation is ^ ~ ^ <S 1. However, as we will see, T^. produces a significant 
effect, because it makes an important contribution to the orientation of the magnetic axis. 

Returning now to the angular variables, we can use a first order approximation in a 
and consider the quasi-stationary regime. To simplify the expressions, instead of dealing 
with the very similar variables Q and /9, we will replace the last angle by a new angle 
7 defined as in Fig. 1. It satisfies cos/3 = cos a cos ^ + sin a sin ^ cos 7, and we have 
Oo-f2eXM = sin 7 sin a sin 6'. Thus, with this substitution and at first order in a, assuming 
^ ~ 1, — = we have: 

Q — — ujf a cos 'f — ujyz sin 6 cos 6 — — f)Q;sin7, (18) 

-'c 

a — —cUxCOsOsinO sin^ — cufa — cUyzCOS^sinOcosO , (19) 

7 = —Q. ( 1 — ^ ) — — cos 7 sin 6 cos 6 + sin 6 cos 6 sin 7 . (20) 
V IcJ Oi a 

The solutions for a and 7 can be decomposed in a transient dynamics, with a characteristic 
time T/, plus a slow varying time-function. This implies that a and 7 will reach a 
quasi-stationary regime in a few seconds. Prom here on they will be driven by the slow 
time-dependence of 6. This assertion can be verified by evaluating the Liapunov exponents 
at the equilibrium point whose real parts are —oJf. The equation for 9 contains only small 
frequencies. The explicit solutions at the fixed point for a and 7 are: 

sm7 = ± \ / ^ ^21) 



a 



"^-^"^y^ sin ^ cos (22) 



These results arc consistent with our previous discussion. In particular a becomes of 
order ^^'^j^ , and thus the approximation a ^ 1 is totally justified. We can also see here 
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that the gravitational dragging torque, despite its noticeable magnitude according to eq. 
(S), has the only effect of renormalizing the angular velocity f2 by a correction of the order 



of 10%. In what follows we use k = 1 — 



0.9. 



Substituting the expressions (|2T|) and (^) into eq. (|T8[), we finally obtain 



6^x00 fC 

3 Ic c3 ujj + fi:2^]2 V 5Qm 
and we have for the angular velocity of the star 



KUyz I sin 9 cos 9 



2M' 



IVl / 



Vyz sin^ 9 + i)yz]^^ ■ 



yz 



With these results, the expression for the breaking index becomes: 



n 



!J2 



QuxUJfC 



KUyz ] sin 9 cos 9 (^Uy^ sin^ 9 + Vy^ 



d9 



(23) 



(24) 



(25) 



The equations of motion (^), (|2^) , and the expression (|25|) for the breaking index show 
that the effects of the magnetosphere are indeed relevant to explain the evolution of the 
magnetic dipole and the angular velocity of a neutron star. 



4. Discussion 



The model depends on several parameters. The mechanical characterization of the 
system is given by h/Ic and R. The lower magnetosphere is described by Vx-, whereas the 
upper one is represented by Vyz and Vyz. The effective friction between the core and the 
crust is given by uof. We have observational information on fi, VL and n. Furthermore, we 
expect to have a star radius of the order of i? ~ 10 km, a ratio of the total momentum of 
inertia and the crust momentum of inertia h/Ic ~ 10^ and k ~ 1. The angular velocities of 
neutron stars are in the range 1 s^^ > Vt > lO'^s^^, whereas a reasonable value for uof is of 
the order of ls~^, and thus we can assume that p^2o2 ~ 1- 
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The simplest situation we can consider corresponds to a star without a magnetosphere. 
In this case = ^yz = ^ and = 0, and hence we have: 



n = -;;^sin^^, (26) 



n 



(28) 



The first equation states that Cl is always negative, which is consistent with the fact 
that the star is losing energy by electromagnetic radiation, and thus the angular velocity 
is constantly decreasing. The second equation implies that the magnetic dipole slides 
toward the direction of the axis of rotation if < 1 or to the equator if > 1, and 
hence the breaking index, as is shown by the last equation, becomes greater or smaller 
than 3 respectively. For example, if we consider the Crab pulsar, for which Q = 190 s~^, 
h = -2.4 lO-^s-^and n = 2.5, they lead to: 

^i^^SlO^sin^^t/r, (29) 
a;/ = 1.2 (1 + 1.7 10"^ tan^ O) s'^ , (30) 
^~ 10"^ tan^yr-^ (31) 



The resulting value for Uf, of the order of seconds, is in good agreement with the 
theoretical expectations and the bounds derived from the Vela and the Crab pulsars. But 
this value for cuf in fact corresponds to a fine tuning to have — 1 ^ 10~^, which is 
rapidly spoiled as Q decreaces. The effects of the electromagnetic aligning torques are much 
higher than the corresponding ones for a rigid star. In particular for a young star they 
give a characteristic time of the order of 10 yr. For this reason the magnetic moment of a 
neutron star in the vacuum will reach the rotation axis or the equator, and stabilize in a 
rather short time compared with the age of the star. If it tends to the rotation axis the 
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breaking index grows much bigger than 3, and if it falls in the equatorial plane the breaking 
index becomes exactly 3. Clearly, this is not the situation observed in the known pulsars. 

The effects of the magnetosphere introduce a qualitative change in the dipole alignment 
behavior. The lower magnetosphere can be modeled by the Goldreich- Julian currents, in 
which case we can take = —1/4. At the moment there is no suitable calculation for the 
non anomalous torques, generated by the upper magnetosphere, and thus we will maintain 
the corresponding coefficients as phenomenological parameters. Then we have: 

n = ----^(uy,sm''e + uy,)n^ , (33) 



n — 3 



lit f SuJfC \ . n n d / .On ~ 

' - 3l llOQ^ + '^"^ M ^ + 



(34) 



There is a very important remark to be made about the Q behavior. Taking into 
account that we have —^j^ — the evolution of Q when -^^^ + z/j^^ ~ 1 is dominated by 
tan 6* oc e~f^ . For a typical neutron star ^ ~ 10~^rn, and thus in a very short time the 
magnetic dipole lines up with the rotation axis, and from there on ^ = 0. This behavior is 
a consequence of the relative freedom of the crust respect to the core, which increases the 
velocity of alignement by a factor But this is not the only possible behavior. Another 
one can be realized if there is an equilibrium point for the dynamics of Q at 

^.(«) = . (36) 

satifying v'yj^ff) cos 6* > to be stable, where the prime indicates a derivative with respect 
to Q. In this case Q will rapidly adjust to the equilibrium value and its dynamics will be tied 
to the angular velocity dynamics 

i>—^K^A%\i^, (36) 
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which implies that 9 slides towards the rotation axis. Besides this, from eq. (35) we have 



If Pyz ^ 100,2^ sin^ 6 ~ —T-'yz sin^ 6', the angular velocity VL decreases throughout the life 
of the star; otherwise, if i)yz < Y^^sin^6', the surroundings accelerates the star. The first 
situation seems to apply to the known pulsars. 



There are theoretical arguments suggesting that ujj depends on VL^, with k < 2 (|Alpar 



Langer fc Sauls 1984| ; [Mendell 1997|) . The angular function VyziO) is of order one, but during 



the life of the star the angular velocity constantly decreases. Thus, at a given moment the 
equilibrium point condition (|35|) cannot be maintained any more. When -^^^ becomes 
significantly greater than Vyzid) there is a change of regime and the dipolar moment rapidly 
aligns with the rotation axis. For a young star evolving at the equilibrium point, where Vy^ 
is of order one, we can estimate a lower bound for tj because it must be close to or greater 
than Yo^^- For example, for the Crab pulsar it is Tf ^ 0.25s and for the Vela pulsar 

^ 1.8s, consistent with the known upper bound. In this regime the breaking index is 
given by 

77 (i^vz sin^ + i^'uz) 

n-?, = {k-2)'-^ y • 2. { - (38) 
^yz [I'yz sm + i'yzj 

which has the correct order of magnitude. To have n < 3 it must be 
(^yz sin^ 9 + Vy^ cos6 < at the equilibrium point. 

Although the aim of this model is to describe the overall evolution of the magnetic 
dipole, a question that naturally arises is what it can say about the glitches, which up to 
this point have not been considered. If we suppose that when a glitch happens there is a 
change in 6, from eq. (|33|) we have 

An [vyz sin^ + Vy:,) 



{uyz sin^ 6 + Pyz 



■AO . (39) 
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The magnetosphere-dependent factor can be assumed of the order of the unity, whereas 
the change in the angle 9 is at most of the same order as j^-^, where AQc is 

the crust angular velocity jump just after the glitch. Thus this relation implies that the 
relative change of (l is at most ^ ~ j*-^. For example, in the case of the Crab pulsar 
this relation tells us that ^ ~ 10^"^ — 10~^, whereas for the observed ghtches it is 
10^*^ ^ ^ ^ 10^'^. This shows that these phenomena involve some non systematic and 
probably very complicated factors, such as deformations, earthquakes and breakings of 
the crust, which are not considered in the present model. For this reason it gives only a 
lower boundary for the ratio between the increment of the angular acceleration and the 
acceleration itself. The change of the 6 and ujf parameters during the glitches could produce 
a departure from the equilibrium point that would be reached again after a transient time 
of the order of ten years. 

5. Concluding remcirks 

In this paper we have developed a two-component model for the evolution of a neutron 
star in presence of a magnetosphere. It considers the core and the crust with a fixed 
magnetic dipole moment, taking into account the normal and anomalous torques that act 
on the dipole and the friction between the core and the crust. We also anahze the effects 
of gravitational corrections and show that they only introduce a renormalization in some 
parameters, but do not affect the qualitative behavior of the system. We have solved the 
equations of motion taken advantage of the very different characteristic time scales that 
emerge from the complete dynamics. 

The anomalous torques are usually supposed to be irrelevant for the alignment of the 
magnetic dipole, but here we show that in fact they are very relevant. This is because the 
crust-core interaction allows the angular momentum of the crust to evolve independently 
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of the magnetic moment by interchanging angular momentum with the core. Despite this 
quahtative change in the ahgnment, the energy loss due to dissipation is very small and the 
equation governing the angular velocity is the same as for a rigid star. 

Another effect of the crust-core interaction is an amplification of the classical alignment 
velocities by a factor that greatly reduces the alignment times. This result and the 
observations lead us to the conclusion that these torques are not directly governing the 
dynamics, but there is an equilibrium point which effectively drives the angle 6. This 
equilibrium is reached by the interplay of the aligning effect of the x axis torque, which 
depends on the friction between core and crust, and the y — z plane aligning torques, which 
depends in a crucial way on the magnetosphere. The equilibrium point naturally moves 
with the caracteristic times of the dynamics of the angular velocity and leads to a breaking 
index near 3, as is actually observed. This equilibrium point regime works for rapidly 
rotating young stars. When the angular velocity becomes small it cannot be established 
and the magnetic moment will rapidly reach its final state. The equilibrium point regime 
requires a crust-core friction in reasonable agreement with the theoretical expectations and 
the observational boundaries. 

One can be tempted to extend this model to the study of the glitches, but the results 
do not agree well with the observational data from the Crab pulsar. We obtain only a 
lower boundary for This is to be expected, because in these phenomena it is very likely 
that more complicated processes, such as crust deformation, fractures and earthquakes, not 
taken into account here, have a very significant role. This suggests that other effects could 
be included in this framework. For the electromagnetic interaction, the actual effect of 
anomalous torques due to quadrupolar magnetic moments should be investigated, although 
their average effect is expected to be strongly suppressed in a rigid star model ([Good fc JNg 



1985). The effects due to crust elasticity and oblateness could also be important as analized 
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in references ( Goldreich 1970 ; |Macy 19"73 ), and they might be relevant to extend this model 



for describing the glitches dynamics. 

In summary, this model gives a consistent description of the overall evolution of the 
main parameters of a neutron star. As long as the magnetic moment can be considered 
constant it is a suitable aproach for understanding the magnetic field and angular velocity 
beheavior and, as has already been said, it can be considered a first step for constructing 
a reliable description of the global dynamics of a neutron star. In particular it reconciles 
the fact that the angular velocity can be permanently decreasing by radiation, and that 
simultaneously for a young and rapidly rotating star we have a breaking index n ^ 3. For 
an older one the magnetic dipole aligns with the rotation axis and stabilize. 
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Fig. 1. — The crust angular velocity fic defines the z-axis. The magnetic moment M is in 
y-z plane, with an angle 6 with respect to de z-direction. The core angular velocity JIq forms 
an angle a with Oc and {3 with M. 
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